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We calculate the second and third virial coefficients of the effective sphere-sphere interaction due 
to polymer-induced depletion forces. By utilizing the anisotropy of a typical polymer conformation, 
we can consider polymers that are roughly the same size as the spheres. We argue that recent 
experiments are laboratory evidence for polymer shape anisotropy. 
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Over sixty years ago, Kuhn studied the conformations 
of polymer chains |lj and recognized that typical con- 
formations of ideal chains are not spherically symmetric. 
The intuitive idea of a symmetric shape is a result of the 
isotropic end-to-end vector distribution of a random walk 
spherical symmetry results from implicit rotational 
averaging of the polymer. In fact, a typical polymer con- 
formation is anisotropic, with an aspect ratio of roughly 
3.4: 1 |. 

Nonetheless, there is little laboratory evidence for this 
asphericity. In solution, polymers rotate randomly and 
there are no significant polymer-polymer correlations. 
Thus small angle light scattering measurements of poly- 
mer solutions can only determine the average princi- 
pal axis of the polymer shape. If there were a way to 
induce strong orientational correlations, light-scattering 
measurements could discern an asymmetry. We will 
show that the induced attraction between regions of de- 
pleted polymer concentration (inclusions) is a probe of 
this shape anisotropy. 

The shape tensor characterizes the spatial distribution 
of monomers: 
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where R a (n) is the position of monomer n, a and (5 label 
the Cartesian coordinates and R a = jj J Q dnR a (n) is 
the polymer center of mass. The polymer radius of gy- 
ration is simply Rq = Tr(M). Moreover, the eigenvalues 
of M a j3, \\ < \\ < A3, are the average squared radii of 
gyration along the principal axes of inertia. Simulations 
HQ have determined that the the most likely shape has 
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Indeed, the shape asymmetry exists and is rather large 
1^1 . Exploiting this shape anisotropy as a calculational 
tool is the main theoretical component of this letter. 

Though hard spheres only interact by direct contact, 
entropic effects of other particles present can induce long- 
range interactions. These sorts of forces are responsible 
for liquid-crystalline order in lyotropic systems Q and 
surface crystallization in hard-sphere fluids [Q. It is in- 
structive to consider the virial expansion for a gas of 
identical balls. Around each sphere of radius r there is a 
sphere of radius 2r from which the centers of the other 



spheres are excluded. When two spheres are close their 
excluded regions overlap leaving more free volume for the 
remaining spheres and hence a larger entropy. Thus the 
propensity for "close" configurations can be interpreted 
as an entropic "depiction" force. 

Recently ||, monodisperse polymers (specifically, A- 
phagc DNA) have been used to induce depletion forces 
between polystyrene spheres. To model this system at 
low concentrations, one might replace the polymers with 
spheres of radius Rq, the polymer radius of gyration. 
Asakura and Oosawa derived a simple formula for the 
potential between two large spheres of diameter a in a gas 
of smaller spheres of diameter D. The Asakura-Oosawa 
(AO) potential is S: 
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where R is the distance between the centers of the large 
spheres, v is the volume of a small sphere, II is the os- 
motic pressure of the small-sphere gas and A = 1 + D/a. 
The above approximation simply sets D = 2Rq- This is 
obviously a crude approximation to the true system. A 
complete analysis at length-scales longer than the poly- 
mer persistence length (50nm) would count the num- 
ber of self-avoiding random walks which avoid the two 
polystyrene spheres. 

In general, the effective potential is of the form [/ c ff = 
IiV(R) where II is the osmotic pressure and V is an 
R dependent recovered volume. At low concentrations 
the osmotic pressure is not adjustable: II = k B Tc. The 
physics all lies in V(R). The AO model gives a one- 
parameter family of functions, depending on the effective 
hard-sphere diameter D. In principle one could derive a 
virial expansion for this potential with each term involv- 
ing the evaluation of a set of cluster integrals, each of 
which involves integrations over polymer degrees of free- 
dom. We will derive a different one-parameter family 
based on the known shape distribution of polymers and 
argue that the data in || is the first laboratory evidence 
of the conformational anisotropy. 

We start by calculating the classical configurational 
integral Qn(R), the sum of Boltzmann weights over all 
conformations of N polymers with two spherical inclu- 
sions separated by R. The sphere-sphere effective poten- 
tial U c ff(R) is: 
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Although Qn(R) has terms which are independent of R, 
the resulting effective potential U e s(R) does not: they 
cancel between numerator and denominator in (Q). 

To calculate Qn(R), we sum over all the conformations 
and placements of N polymers with the excluded- volume 
Boltzmann weights: 1 if the polymers and the spheres do 
not overlap and otherwise. We split the integration over 
each polymer into three parts. The first is an integration 
over the center of mass, the second an integration over 
all rigid rotations of the conformation, while the third is 
the remaining integration over "internal" degrees of free- 
dom. This final integral is over each unique conforma- 
tion - two conformations are equivalent if one is merely 
a rigid rotation or translations of the other. We will only 
integrate over one representative from each equivalency 
class. Denoting the space of all such "internal" polymer 
conformations by T, we have: 
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where G?r^i nt is the measure on the space of conforma- 
tions for polymer i, r< is its center of mass and fij is 
its rigid rotation. We further divide the integration over 
T by characterizing each polymer conformation by its 
principal axes. Defining g(Ai,A2,A3) is the number of 
conformations with axes A, , we have 
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To pass from (||) to @ we assumed that the internal de- 
grees of freedom did not affect the interaction potential 
U. This is, of course, not precisely correct. Our approx- 
imation replaces each polymer by a solid ellipsoid, and 
then considers the potential due to only to this shape. 
While this certainly removes many degrees of freedom, it 
includes more degrees of freedom than replacing the poly- 
mers by spheres. The proof of the pudding shall be in 
the eating - we will see that this approximation is valid 
by comparison with data Thus our approximation 
replaces the monomer-sphere and monomer-monomer 
potential with a sum of pairwise, ellipsoid-ellipsoid or 
sphere-ellipsoid terms. Each term is infinite for any over- 
lap and zero otherwise. 

We now reduce the complexity of the integration in 
@ by using the shape distribution of polymers. Since 
the distribution of polymer shapes is peaked around 
the prolate spheroid, we only consider those polymer 
shapes. This approximation does not account for the 
entire space of principal axes, though we believe that 
it does characterize the polymer conformations better 



than a sphere. More importantly, our approximation al- 
lows us to consider polymers which are roughly the same 
size as the included spheres. This has a great advan- 
tage when comparing to the experiments of Verma, et al 
M. By comparison, most work in this field has treated 
the inclusions approximately while correctly modelling 
the medium as a gas of random walkers. In |l0| field 
theoretic methods were used to obtain the depletion po- 
tential. There, in order to calculate reliably, the au- 
thors considered the interaction between inclusions much 
smaller than the polymers, i.e. a <C Rq- in the 
opposite extreme, one might consider polymers which 
are much smaller than the spheres. In this case, it 
is appropriate to study the induced Casimir force be- 
tween two walls and, to then approximate the interaction 
between spheres via the Derjaguin approximation 
By approximating the polymers as prolate spheroids, 
we can treat the inclusions exactly - a great advan- 
tage when the two extreme limits are not applicable. 
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FIG. 1. Sequence of approximations. The actual polymer 
performs a self- avoiding random walk (a) which has a typical 
prolate spheroidal shape (b). For computational simplicity, 
we neglect the anisotropy between the two shorter principal 
axes and build up the resulting ellipsoid out of overlapping 
spheres (c). 

We have thus reduced the configurational integral Q n 
over all polymer modes to an integral over the allowed 
locations and orientations of a prolate spheroid. In order 
to reduce the phase space somewhat and reduce our com- 
putation time, we will assume that Ai = A2. In the case 
of interest, we consider a prolate spheroid with aspect 
ratio Vll-8 : 1 ~ 3.4 : 1.0, so that 
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where Rq is the radius of gyration of the polymer. To 
facilitate the numerical evaluation of integrals, we con- 
struct the ellipsoid out of overlapping spheres. Figure 
1 depicts our sequence of approximations. We note that 
choosing A3 S> A 2 = Ai would reduce our analysis to that 
in which considered the depletion force between two 
spheres induced by thin rods. 
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volumc potential between particles i and j, and labelling 
the spheres A and B, Qn(R) can be rewritten as 
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The product in ^| is a sum of terms which can be grouped 
by the number of polymer positions and rotations that 
are freely integrated over. The first term is proportional 
to (4ttV) n , the configurational integral for free ellipsoids, 
where V is the volume of space minus the volume of the 
two included spheres. Subsequent terms have fewer pow- 
ers of V. Since the polymers are identical, these correc- 
tions include combinatoric factors involving N. We take 
the limit N — * 00 and V — ► 00 keeping c = N/V constant 
to find the virial expansion in c. 

It is both convenient and instructive to graphically rep- 
resent these "cluster" integrals |13 in terms of Mayer 
cluster graphs. The first two terms in U e s(R) are: 
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where the open dots represent the spherical inclusions 
and the closed dots represent the ellipsoids. The inte- 
grals in (||) are difficult to compute analytically and thus 
were evaluated numerically via a Monte Carlo algorithm: 
10 3 different angles and 10 4 different points were chosen 
in a volume which included both spheres and which did 
not exclude any possible orientation or location of the 
ellipsoids. We calculated these integrals by this random 
sampling weighted by the appropriate phase-space vol- 
ume factor. 

To compare with experiment ||, we took the sphere 
diameter to be D = 1.2/xm and Rq — 0.5/im. Know- 
ing Rq enables us to find the length of the ellipsoid, 
L = 3.4 x 2R G /Vl3£ w 0.92/zm. We have chosen the 
hard-ellipsoid radii to be the mean square radii of gy- 
ration, which is possibly naive: the hard-ellipsoid size 
should only be proportional to Rq. Since, in a random 
walk the density decays as 1 jr there is no natural length 
scale which cuts off the excluded- volume interaction. In- 
deed, light scattering experiments |yj have found that 
the effective hard sphere radius is roughly half the ra- 
dius of gyration. The relation between Rq and the hard- 
ellipsoid size must be determined through the depletion- 
force experiments we are modeling. 

In Figure 2 we plot our results as a function of con- 
centration. Until one considers concentrations near the 
polymer overlap concentration c* = l/(4irR G /3) w 
2//im 3 the third virial coefficient, responsible for a 
repulsive "anti-correlation hole" , is a small perturba- 
tion to the leading term in (^). Thus if we re- 
strict our study to the dilute polymer regime, the 
leading term in the virial expansion is sufficient. 
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FIG. 2. Numerical result for the effective potential U e a for 
differing concentrations and Rg = 0.5/xm. We have kept 
terms up to second order in the polymer concentration. Note 
that at the highest concentration there is a small repulsive 
bump at R « 2(im. However, for c < c* ~ 2//xm 3 , the third 
virial coefficient is a small contribution to U e g. 

We compare our model with the data Q and the AO 
model at c = 0.5//im 3 for two reasons. This concen- 
tration is below the overlap concentration c* and at the 
same time is large enough that the well depth is on the 
order of k B T so that data can be reliably obtained. We 
can adjust the effective radius Rg for both the AO model 
and our ellipsoid-based model. We find that a good fit 
results for the AO model with Rq — 0.42/im and for our 
model with Rq = 0.8/xm, in comparison with the light- 
scattering-obtained value Rq — 0.5/im. In Figure 3 we 
show the data along with these two one-parameter fits. 
We have checked other concentrations and have found 
that at c — 1.0//im 3 the theory and experiment also 
compare favorably with the same effective radii. What 
should one conclude from this agreement? It is possible 
to interpret the data as arising from either shape simply 
by adjusting the size of the shape. However, one should 
start from a microscopic picture based on the polymer 
physics of the allowed chain conformations. Only from 
this perspective can one properly interpret the data. 

Since the radius of gyration of the A-DNA can be cal- 
culated from its molecular weight and persistence length 
to be Rq = 0.5/im we can consider two different zero 
parameter fits: our model and the AO model. The re- 
sult is shown in Figure 4. Note that the data lie between 
our calculation and the AO model, suggesting that one 
could smoothly deform the AO sphere into an ellipsoid 
and find a best fit aspect ratio to fit the data, giving 
a one-parameter fit. Indeed, we have made a number 
of limited runs on the fully anisotropic shape satisfy- 
ing (||). This results in a curve which is roughly 30% 
deeper than that in Figure 4, which is better than the 
simple sphero-cylinder result. Finally, we have made 
similar comparisons between theory and experiment at 
c = 0.1/nm 3 , c = 0.2/^m 3 and c = 1.0/fim 3 . We find 
that at c = 1.0//im 3 the comparison is similar to that 
shown in Figure 4. At the two lower concentrations, 
where the data is difficult to collect, neither our model 
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nor the AO model make very good predictions. 
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mations by treating them as rigid ellipsoids is appropriate 
for static properties, it is not clear at all that the dynam- 
ics reflects this. In particular one might ask whether a 
polymer ellipsoid rotates to a new orientation slower or 
faster than it deforms into that orientation. Finally, our 
analysis could also be used to study the depletion inter- 
action by actual ellipsoidal objects, such as bacteria Jl5| . 
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FIG. 3. Our model and the AO model are fit to the exper- 
imental data for c = 0.5//im 3 . We varied only the effective 
Ra as a parameter. The fit gives Rg = 0.42 and 0.8 for the 
AO model and our model respectively. 
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FIG. 4. Comparison of the AO model, the model presented 
here and the data at c = 0.5//im 3 . For both our model and the 
AO model we have taken the theoretical value of Rg ~ 0.5pm. 
There are no free parameters in the models. We also plot the 
average of the AO model and the ellipsoid calculation. 

Returning to the original question, are polymer shapes 
anisotropic? Though the full distribution g(Ai, A2, A3) is 
peaked, it has a finite width. We can incorporate this 
width by choosing an appropriately weighted admixture 
of shapes. Note that we could easily get a very good fit 
just by mixing spheres and ellipsoids - the mixed second 
virial coefficient is just a linear combination of the two 
"pure" coefficients. It is easy to see from Figure 4 that 
a 50-50 admixture of spheres and ellipsoids would give a 
remarkably good fit to the data, without any adjustable 
parameters. While we have no basis for this weighting 
of shapes, we nonetheless take this as direct evidence for 
the anisotropy of polymer conformations. 

In closing, we note that while counting polymer confer- 
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